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The spherical field formalism — a nonperturbative approach to quantum field 
theory — was recently introduced and applied to 4 theory in two dimen- 
sions. The spherical field method reduces a quantum field theory to a finite- 
dimensional quantum mechanical system by expanding field configurations in 
terms of spherical partial wave modes. We extend the formalism to 4 the- 
ory in three dimensions and demonstrate the application of the method by 
analyzing the phase structure of this theory. 



I. INTRODUCTION 

Recently a new nonperturbative approach to quantum field theory, called spherical field 
theory (SFT), has been introduced by Lee [0. The idea of this approach is to reduce a 
quantum field theory to a finite-dimensional quantum mechanical system by expanding field 
configurations using spherical partial wave modes. The dimension of the system is made 
finite by truncating the partial wave expansion by neglecting high spin waves. The resulting 
quantum mechanical system can then be analyzed numerically in Euclidean spacetime with 
the use of diffusion Monte Carlo methods (DMC). 

SFT was developed in the context of (0 4 )2 theory in [[!]]. The purpose of the present work 
is to demonstrate that SFT can be extended beyond two-dimensional theories by developing 
the formalism for (0 4 )3 theory. In fact, we start out with a general extension of |TJ to an 
integer number of dimensions N > 2 and specialize to N = 3 when we renormalize the 
theory. 

The paper is organized as follows: We begin by expanding the field configurations in 
the generating functional for 4 theory in N dimensions in terms of spherical partial waves 
and integrating over the angular coordinates in the resulting expression. Then, allowing the 
radial coordinate to play the role of a time coordinate, we use the Feynman-Kac formula 
to write the generating functional in terms of a Schrodinger time evolution operator called 
the spherical-field hamiltonian. Next we apply functional differentiation to the generating 
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functional to yield the free two-point correlators of the partial wave fields. Then we discuss 
the renormalization of (0 4 )3 in the context of spherical field theory. Finally, to illustrate the 
application of this new method we model the quantum mechanical system defined by the 
spherical-field hamiltonian using DMC in order to analyze the phase structure of (0 4 )3- 

Consistent with Magruder's study of (</> 4 )3 theory || we find that the theory undergoes 
a transition from a broken-symmetry phase to a symmetric phase if the mass-squared pa- 
rameter /i 2 is negative. In the case of fi 2 positive Magruder's analysis does not determine 
if spontaneous symmetry breaking occurs in this theory. Several authors have investigated 
this question using various nonperturbative methods; the results of their analyses are, how- 
ever, in disagreement. Cea and Tedesco || studied a generalization of the Gaussian effective 
potential in which they calculated the two-loop corrections to the Gaussian approximation. 
They found no spontaneous symmetry breaking. In contrast, Stancu |4j] used an expansion 
of the effective potential based on the variational Gaussian effective potential method and 
found that a phase transition does exist. Peter, Hauser, Thoma and Cassing |J also found 
spontaneous symmetry breaking by applying the coupled set of equations of motion for 
Green functions up to the 4-point level truncated by the use of the cluster expansions for 
n-point Green functions. With our analysis we find that the theory does not undergo a phase 
transition if fi 2 is positive — the vacuum remains invariant under the reflection <p — > —cf). 

For the phase transition observed for fi 2 < we compute the critical coupling and the 
critical exponents v and (3. Recent calculations of the critical exponents for the three- 
dimensional Ising universality class — to which (0 4 )3 is believed to belong — have been per- 
formed to a very high level of precision using various approaches such as Monte Carlo sim- 
ulations of 4 lattice models (MC) ||, high temperature series expansions (HT) J7], d = 3 
expansions and e-expansions ||. At the present early stage in the development of SFT we 
are limited by the underdeveloped state of the Monte Carlo algorithms available to model 
the dynamics of systems governed by spherical-field hamiltonians and so we do not intend 
our estimates of v and (3 to be competitive with the state-of-the-art results. Rather, we 
compare our estimates with them in an attempt to build confidence in this new method. 

There are important advantages to the spherical field formalism. The method maintains 
exact rotational invariance and it can be generalized to other geometries as was recently done 



with a study of quantum field theory on the noncommutative plane ||10|| . Most importantly, 



it is a continuum approach rather than a discrete approximation method and so it is free of 



certain problems which arise in lattice calculations such as fermion doubling |TT 



2 This state of affairs is not restricted to SFT. For example, see || for a relevant discussion in the 
context of hamiltonian lattice gauge theory. 
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II. THE GENERATING FUNCTIONAL 



The first step in the spherical field theory approach is to perform a spherical partial wave 
expansion of the fields in the generating functional.^ In iV-dimensional Euclidean space the 
generating functional for <p A theory is given by 

Z[J] = Jv4>exp || d N x [C + J4>]} (1) 

where 

i N 2 

£ = -~£ W) 2 - ^ 2 - f/. (2) 

We express the cartesian coordinates in terms of hyperspherical coordinates, 

Xi = t sin t?! sin $ 2 • • • sin i? ^-2 cos (p 
X2 = t sin i?i sin $2 • ■ ■ sin i?j\r-2 sin 93 

X3 = i sin^i sin $2 • • • cos #^-2 (3) 



£jv-i = t sin 1?i COS i? 2 
Xat = tcos^i, 

and decompose the field and the source J into partial waves: 

= E ^>A,J^A,J^ ( 4 ) 

A.m 



A,m 

We have written A for the collection of indices Ai, A2/ ■ ■, A^_2 and Q for the angular 
variables $1, $2, • • • , $n-2, f- Here, and throughout this paper, the sum over A runs over all 
integer values of the Aj such that 

Ai > A 2 > A 3 > ••• > \ N -2 > (6) 

and the sum over m runs over all integer values satisfying 

- Ajv_2 <m < Aiv_2- (7) 



3 It should be noted that covariant Euclidean quantization, an important component of the spher- 



ical field approach, was introduced in |12[. 
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The m are the iV-dimensional hyperspherical harmonics^] which reduce to the usual spher- 
ical harmonics for N = 3. For a useful reference on hyperspherical harmonics as well as the 
hyperspherical Bessel functions we will need below, see [|13[ . 

To express the generating functional in terms of the partial waves we make a change of 
integration variables in the functional integral of Eq. ([[]) from <fi to the set {4>\ m }- The 
Jacobian of this transformation is 



det 



5<f){x') 



det 



(8) 



which is independent of the (f)\ (t). Thus this change of integration variables simply intro- 
duces an overall constant which will be irrelevant for our purposes. 

Now, having inserted Eqs. and (|^) into Eq. ([!]), we can simplify the expression for 
Z[J] by using the orthogonality condition of the hyperspherical harmonics, 



(9) 



to integrate over Q in the action. Before we can do this we must eliminate the angular 
derivatives of the m which appear in the kinetic term of C. This can be accomplished by 
integrating this kinetic term by parts to yield 



1 r N 

^ / d N x £(a, 
A J i=i 



- I d N x 



(10) 



where A denotes the generalized Laplacian operator which can be written in terms of the 
generalized angular momentum operator A 2 as 



A 



i d 

t^dt 



N-l 



dt t 2 



(ii) 



Now the angular derivatives have been removed since the hyperspherical harmonics are 
eigenfunctions of A 2 satisfying the relation 



A2Y \JV)=l(l + N-2)Y x jn) 



(12) 



where, for ease of notation, we have written t for Ai, the first component of A. 

In performing the integration over Q of the interaction term we encounter an integral of 
the product of four hyperspherical harmonics. Let us make the definition 



Ai A2 A3 A4 

mi m'2 m:j m,4 



n ~ T(N/2) 



dVt Yx (9)Yx (9)Yx (n)Y\ (SI) (13) 



where the prefactor to the integral, which is the area of an A^-dimensional unit sphere, is 
introduced for later convenience. This integral can be calculated by using the expansion 



4 We use the phase convention YZ (O) = (-l) m Y^ _ m (^). 
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^Ai,mi (^)^A 2 ,m2 ^> 



T(N/2) 

27rAf/2 C N (\ l ,m l ;\2,m 2 ;\,m)Y Xm (a) 

A,m 



(14) 



where again a prefactor has been introduced for convenience. The coefficients, 
Cn(\i, mi, A 2 , m 2 ; A, m), have been worked out in |TJ]]. They vanish unless m± + m 2 = m 



and so it follows that 
for iV = 3 we have 



Al A2 A3 -^4 

mi ni2 1713 m.4, 



N 



vanishes unless mi +m 2 +7713+7714 = 0. In particular, 



h i 2 h h 

mi ni2 7713 m4 



( — \ \mi+m2+m3 r 

V 1 J u mn,—mi 



-m 2 -m 3 C 3 (£ 1 ,m 1 ;£ 2 ,m 2 ;£,m 1 +m 2 ) (15) 

e=\h-h\ 



xC 3 (£, mi + m 2 ; £ 3 , m 3 ; £ 4 , m 1 + m 2 + m 3 ) 



with 



C 3 (£ 1 ,m 1 ;£ 2 ,m 2 ;£,m) 



(24 + 1) (2£ 2 + i; 



M (2* + 1) 

x(4,4;0,oK 1 ^a;^o) 



(£i,£ 2 ;mi,m 2 |£i,£ 2 ;£,m) 



(16) 



With these results we are in a position to write down our final expression for the gener- 
ating functional. We have 



Z[J] oc 



n^A )m J^p{/°°^4 ph } 



(17) 



where £ j?* 1 is the spherical-field lagrangian, 



,S P h 

-J 



+ 



eh: 

A,m 

flH 2 + i 



^- ^A-m^A, 
2 dt dt 

' N — 2) „_3 



(18) 



^ < ^\,—m < fi\,m ^ *-^A,— m^A 



_„j-N-i y^ y^ y^ y^ I" A x A 2 A 3 A 4 

1 / ; / ; / ^ / ; [ mi ni2 rrt3 — mi— 013 

Ai,mi A 2 ,m2 As,m3 A4 

^ Ai ,mi ^ A2 ,m2 t A3 ,m3 ^ Az 



N 



V4,— mi— m.2— m.3 



and the constant c is defined as 



c = 



9 T(N/2) 
4! 27r"/ 2 



(19) 
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III. THE SPHERICAL-FIELD HAMILTON! AN 



With the Feynman-Kac formula we can interpret the functional integral Eq. (|T7|) as a 
time evolution equation where the radial coordinate t now plays the role of a time parameter. 

The Feynman-Kac formula, generalized to include the class of time-dependent inter- 
actions relevant here, connects the path integral representation of Euclideanized quantum 
mechanics to the Schrodinger representation, 



<t>{t F )=x F 



where 



dt 



ti 



a(t) ( d<fi 



ft +v '<^> 



oc 



(xp\ T exp |— J dtH(t)^\xi) 

(20) 



H(t) 



d 2 



+ V(q,t). 



2a(t) dq 2 

We may use this formula to write down the Schrodinger representation of Eq. (|17|' 



Z[J] oc 



TeXP {-/ C 



where 



H 



? p V) = E(-i: 



d 



dtH S / h (t)\\a) 







2 ^ 9Q X ,- m 9Q X ,r 



(21) 



(22) 



(23) 



+ 



/i 2 t 2 + 



+ J ^" N _ 3 _ N _ t 

) "A,— m"A,m ^A,— m"A.r 



W^E E EEC 

Ai,mi A2,m 2 A3,m 3 A4 
X ^Ai ,mi ^A 2 ,m 2 ^As ,m 3 ^A4 



A 2 A3 A4 

Wl2 "^3 — mi— m 2 — JTI3 



V4,-rrti— m 2 — m-3 



is a Schrodinger time evolution generator which we will refer to as the spherical-field hamil- 
tonian. Following U we take the quantum state at ti = to be a superposition of all 
possible states with each state equally weighted, 



\ 



| a) = f J] dx X, m \i x X,m>}) 

,m I 



(24) 



In the set {x-^ ,}, the and m' take all integer values satisfying the inequalities (§|) and 
(|7|) respectively. Note that, with this notation, 

(25) 



The proper quantum state for tp —* oo is the ground state, |0). 
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We are faced with the difficulty that the spherical-field hamiltonian is written as an 
infinite series. In practice we cut this series off by neglecting partial waves with spin higher 
than some value we call Jmax- Such an approximation is justified because, as pointed out 
in the high spin partial waves correspond to high tangential momentum modes so that 
a partial-wave cutoff J max corresponds to a momentum cutoff 

A 2 ~ (26) 

Therefore the high spin modes will decouple from the theory, just as other high momentum 
effects, if the theory is renormalized to remove ultraviolet divergences. Below, we will 
describe the proper method of renormalizing (0 4 )3 in the context of spherical field theory, but 
to do so we need to derive an expression for the free two-point correlators of the partial wave 
fields which can then be used to evaluate the spherical-field version of Feynman diagrams. 



IV. CORRELATORS 



Correlation functions for partial waves can be computed from the generating functional. 
We begin by writing the generating functional Eq. ([[]) as a power series, 

w = 1 + \ I ^S^ e " P ' (CE_2/)l7(a;)l7(2/)/(p2) + ' ' • (27) 



'[0] 

where f(p 2 ) is the propagator, 



f{p 2 ) = / d N x jP x <0| 0(*)0(O) |0> . 



(28) 



We insert the partial wave expansion of J(x) given in Eq. (§) and expand the plane waves 
with 



ipx 



(N-2)\\ 



2n N ' 2 
T(N/2) 



(29) 



where k = \p\, t = |cc| and jN, n is the hyperspherical Bessel function of order n, 

T{a)2 a -\J n+a {x) 



JnA x ) 



(N-4)\\t° 



(30) 



where a = Then we can perform the angular integration in Eq. (P7|) to obtain 



Z[0] 



N-l 



(N-2)\\ 
T(N/2) 



E 



dk dt x dt y 



N-1j.N-1j.N-1 



xk iy ~H 



tf- ^(k 2 )J X ,-Jtx)J\Jt y )j N ^kQj N ,e(kt y ) 



(31) 



From this expression we may calculate the two-point correlator for 0^ r 
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Z[0] ViSJxJtJtr^JxJt 



-z[J] 



(32) 



J=o 



2 N-2 



(N — 2)!! 
T(N/2) 



f(k 2 



The free-field two-point 0_\ m correlator is obtained by inserting the free propagator 



^-y into Eq. (|3~2"|). The result is 



m , ,2a 



(JV-2)!!r(a) 



(33) 



2T(N/2) 

x [B(t x - t 2 )k N/ (nti)i Nt i(fit 2 ) + 6(t 2 - t^kuji^iN^inh)] , 

where fc/v,n and are the hyperspherical modified Bessel functions of order n which are 
given by 



k N ,n( X ) 



inA x ) 



T(a)2 a - 1 (N -4)\\t a ' 

n+a ( 



T(a)2 a - l I n+ Jx 



(N — 4)!! t a 



(34) 
(35) 



V. RENORMALIZATION 

Hereafter we will specialize to three dimensions.^ In this case the theory can be made 
finite by adding a mass counterterm to the lagrangian Eq. (^]) which subtracts all primitively 
divergent self-energy diagrams. In three dimensions there are two such diagrams. One of 
these, shown in Fig. |l|(a), is local and so it is straightforward to write down the appropriate 
counterterms corresponding to this diagram in the context of spherical field theory. We 
must explicitly subtract the contribution of all diagrams having the form shown in Fig. |2](a). 
Using Eq. ( ) , and recalling the spherical-field lagrangian Eq. (|T%|), we find that the 
counterterm which corresponds to a diagram of this form is proportional to 



lr'fict 2 



l £' 
-m m' 



ki'(fd)i£>(iit)4>t- m (t)^t, m (t). 



(36) 



The proportionality constant is a product of two factors. One is the degeneracy of the quartic 
term in the lagrangian which gives rise to such a diagram, or equivalently the number of 
distinct permutations of the columns in the bracket symbol. The other is the number of 
possible contractions of this quartic term which correspond to the diagram. 

The other divergent diagram is the so-called sunset diagram shown in Fig. |l|(b). This 
divergence can also be eliminated by using a local mass counterterm. In ordinary field theory 



5 We will no longer write dimension subscripts, i.e. k^ n (x) will be written as k n (x), etc. 



8 



(a) 



FIG. 1. The primitively divergent self-energy diagrams in ordinary perturbation theory. 



:rn 




(a) 



(b) 



FIG. 2. The primitively divergent self-energy diagrams in spherical field theory. 



one finds that the divergence is logarithmic, and so a single subtraction at any value of k 2 
will cancel it. It is natural to choose k 2 = 0. In the context of spherical field theory, we 
work in position space. To obtain a local counterterm we must consider all diagrams of the 
form shown in Fig. |2](b) and we must integrate over the coordinate of one of the vertices, 
say t'. However, we must do this in such a way as to keep our renormalization condition 
translationally and rotationally invariant. 

Working again in ordinary perturbation theory, we consider the sum of the sunset dia- 
gram and its counterterm in position space, 



Here, the first term on the right hand side denotes the diagram itself and the second term 
denotes the counterterm. For convenience, we rewrite the second term as 5a(x, Q x >) ~ 
absorbing the angular part of the delta function into 5a. We want the amplitude to vanish 
at k 2 = 0. To impose this condition in a translationally and rotationally invariant way we 
demand that 




(37) 




(38) 



where f{k 2 ) is some function of k 2 which vanishes as k 2 — > 0. 
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We can move the exponential factor independent of x' to the right hand side of the 
equation and make use of Eq. (p9|) to obtain the condition 

J dV j e (kt')Y ijrn (n x ,)X(x, x') = j,(kt)Y etm (^x)f(k 2 ) (39) 

which must hold for all possible values of I and m. Finally, we take the limit k 2 —>■ 0. In 
this limit, je(kt) oc (kt) and so we have 

J dQ x , Y itm (n x ,)5a(x, n x ,) = i j°° dt't' (i+2) J d{l x , Y e>m (n x ,)a(x, x') (40) 

for all I and m. This relation tells us the proper way to integrate a diagram of the form of 
Fig. |2|(b) over t' in order to obtain the corresponding counterterm. With this information 
and again using Eq. ([33]), we can write down this counterterm. It is proportional to 



1 



■l) m c 2 t 2 



mi ni2 m — m\—m2 —m 



2 



4% 2 AMt)^ m (t)<f>i, m (t) (41) 



where 



4ik/ 3 ( x ) = k h( x ) h 2 (x) h z (x)^j dx' x' {n2) i £l (x')i e . 2 (x')ii> 3 (x') (42) 
^(x^e^x^e^x)^ [ dx' x'^ + ^k ei (x')ke 2 {x')k £3 (x'). 



' ' ~ ' X*- Jx 

Again, the proportionality constant is a product of two factors. Each diagram of the form of 
Fig. ^|(b) arises from the product of a quartic term from the lagrangian Eq. (|18|) evaluated 
at coordinate t and the same quartic term evaluated at coordinate if . Therefore, one factor 
in the proportionality constant is the square of the degeneracy of this quartic term. The 
other factor is the number of possible contractions of the product of these quartic terms 
which correspond to the diagram. 



VI. THE PHASE TRANSITION 

In p|, Magruder proved the existence of a phase transition in (0 4 )3 theory by extending 
Chang's proof |15| for (0 4 )2 theory. Let us briefly recall Magruder's analysis. Consider the 
two lagrangians: 

A- = \{dM^4>) - - W + i^0 2 , (43) 

l_ = \{dM^4>) + -/-<P 2 - W + 

where n 2 + and [i 2 _ are both positive parameters. The mass counterterms, 5[i 2 + and Sfi 2 _, are 
chosen so that the theories are properly renormalized using the same conditions we have 
used above. Expressed in terms of a momentum cutoff A they are defined as 
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1 



. 6 



0.4 



0.2 



0.5 1 
FIG. 3. The functional relation between c 



1 . 5 

c_ 



2 . 5 



and c_ = as given by Eq. ([45|) 



5n 2 + = 12c(A - n + ) + 96c 2 In 1 



5 (j? = 12c(A-/i_) + 96c 2 In 



A 
A 



(44) 



The two lagrangians differ in the sign of their mass terms so that in the weak coupling 
limit the theory defined by C + has a vacuum with manifest symmetry under the reflection 
— > —0 while the theory defined by £_ has a vacuum with broken symmetry. 

The two lagrangians are identical if the parameters of the theories are chosen such that 



- n\ + 5fi 2 . 



-u 2 _ + 5/i 2 



(45) 



The solution to this equation, shown graphically in Fig. |3], gives c + = — as a function of 
c_ = — . This solution has the property that where C + and £_ overlap they are inversely 
related: a strong-coupling theory defined by C_ is dual to a weak-coupling theory defined by 
C + . This implies that the theory defined by L_ undergoes a phase transition from a vacuum 
with broken symmetry to one with manifest symmetry at some finite value of the coupling 
strength, c_ = c CTlt . However, Fig. ^| shows that there is no overlap between the two theories 
for c_ ^ 0.3522. Therefore, if c*i rit is below this value the theory defined by C + will have a 
vacuum with manifest symmetry for any value of c + . On the other hand, if c c ]^ > 0.3522 
then there is a finite range of intermediate coupling strengths for which the theory defined 
by C + will have a vacuum with broken symmetry. Since this approach cannot tell us the 
value of it cannot determine whether or not spontaneous symmetry breaking exists in 
the theory defined by C + . 

We can bring the method of spherical field theory to bear on this question by studying 
the quantum mechanical system defined by the spherical-field hamiltonian. The appropriate 
hamiltonian is found by adding the two counterterms discussed above to the spherical-field 
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hamiltonian given by Eq. (|23]) and setting J = 0. We can do this for the Euclideanized the- 
ory corresponding to £_ as easily as for the one corresponding to C + since the renormalization 
procedure is the same for both theories as indicated by Eq. ([H]). The same counterterms, 
Eqs. (|36]) and ([IT]), apply in both cases. To avoid the complication of working with a 
time-dependent hamiltonian we have chosen to investigate the original three-dimensional 
system on a spherical surface of radius R. Then we may set t — R in the spherical-field 
hamiltonian, rendering it time- independent. The full field theory result is recovered in the 
limit Jmax — ► °o, R — > oo. 

As an example, we display the spherical-field hamiltonian corresponding to C + for J max = 

1: 

1 / d 2 d 2 d 2 d 2 \ fi 2 R 2 2 fi 2 R 2 + 2 



H ~ 2R 2 \dq 2 fi + dq% + dqU + dq^J + 2 % < + 2 + + ^al 

+ci? 2 [g 4 + 6g 2 (g 2 + g 2 sl + g 2 al ) + ~ (g 2 + g 2 sl + g 2 al ) 2 ] (46) 
-6/ici? 2 [k (nR)i (fiR) + 3Ari(/xi2)ii(//i2)] (g 2 + g 2 sl + g 2 al 
+48c 2 J R 2 [(/fioC/iit!) + 9ig,i(A*i2)) <£„ 

+3 teiM) + ;J$,i(A*i2)) feo + 9i 2 , s i + <ai 



where 

qi,si = — y= and q 1M = -j= . (47) 

The spherical-field hamiltonian corresponding to £_ for J max — 1 can be obtained from the 
above expression by letting /i 2 — ► — |/i 2 in the mass terms while leaving the counterterms 
unchanged. Due to their lengths we will not include the corresponding expressions for higher 
values of J m ax but they can be written down in a straightforward way using Eqs. (^3), fl3"B"|) 
and (0). 

Equipped with the spherical-field hamiltonian we can use the method introduced in [TB 



to find the vacuum expectation value and the physical mass as a function of the coupling, 
c. Near the phase transition the behavior of these quantities define the critical exponents (3 
and v respectively. 

We calculate the physical mass in the symmetric phase of the theory by using DMC0 to 
compute the matrix element 

f(t) = (b\q ,oe- tH q ,o\b), (48) 

where \b) is any state that is even under the reflection transformation <ft — > — </>. To see that 
one can extract the physical mass from f(t) insert a complete set of energy eigenstates \i) 
corresponding to energy eigenvalues Ef. 



6 A self-contained introductory presentation of DMC is given in [17]. 
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/(*) = £e- 



-tEi 



9o,o 1 6) 1 1 



(49) 



Now the contribution of the ground state vanishes in this sum because both |0) and \b) are 
even under reflection symmetry while go,o is odd. Therefore, the sum is dominated by the 
one-particle-at-rest state in the limit t — > oo. If we shift the energy scale such that E = 
then for large t 



fit) 



~ e 



-mt 



(50) 



where m is the physical mass. 

Next let us consider the vacuum expectation value. In the broken-symmetry phase of the 
theory there are two degenerate ground states if the size of the system becomes infinite, i.e. 
R — > oo. Let us take these two ground states to be |0 + ), which is non-zero only for go,o > 0, 
and |0~), which is non-zero only for go,o < 0. We choose these such that they are normalized 
to unity and so that under reflection symmetry they transform from one to the other. Now 
our task is to calculate the vacuum expectation value (0 + | <fi |0 + ). Using translational and 
rotational invariance we can write this as 



o- 



cr 



4tt J 



dVL (0 



<j>(t,e,<p) 



cr 



cr 



<?0,0 



0" 



Let us rewrite the operator q in terms of a projection operator, 

/oo 
dxofl xofi \xofi) (xo,o\ ■ 
-oo 



(51) 



(52) 



Now we have 



o- 






o+; 











dxoo Xoo 



47T 



(53) 

since |0 + ) vanishes for 



where the lower limit of integration can be taken to be xo,o 
x ,o < 0. 

The task which remains is to calculate the matrix element which appears in this last 
expression. To do this we need to know the ground state |0 + ). The time evolution operator 
exp (—Ht) operating on any state \b) will project out the state of lowest energy as t — > oo 
provided |6) has a non-zero overlap with this lowest-energy state. Since we are working with 
a system of finite size the vacuum will not be exactly degenerate. For finite R the state that 
is lowest in energy is the symmetric superposition of |0 + ) and |0~), 



°*> = ^(l 0+ ) + l r ))- 



Therefore we have 



|(a; ,o|O s )| 2 = lim 

t— »oo 



(b\e- tH \x 0fi ) (x , \e- tH \b) 
(b\ e~ 2tH \b) 



(54) 



(55) 
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where \b) is chosen such that (0 s | b) is not zero. Now, since |0 + ) = a/2 |0 s ) for £0,0 > 0, the 
vacuum expectation value can be written as 



H 



0- 



2 / dx Q ,o xo,o g(xo,o) 
Jo 



where 



9(xo,o) 



IKo|O s )r 



(56) 



(57) 



We can calculate ^(xo.o) by using DMC to evaluate the right-hand-side of Eq. ([55]) . However, 
we will not calculate (0 + | |0 + ) by integrating 2xo,o5 , (^o,o)- Instead we will take 



0+ 



0" 



x, 



0.0 



(58) 



where 1™ is the non-negative value of xo,o for which g(xo,o) has its maximum value. 
This approach is justified because ^(^0,0) becomes sharply peaked as the size of the system 
becomes large and <?(xo,o) satisfies 



dx ,og(xo,o) 



dx 0y0 g(x ,o) = 1. 



(59) 



This approach converges to the same result as integrating 2x 0i o5 , ( :r o,o) m the limit R — > 00, 



however, as indicated in |L6 it is less susceptible to systematic errors for finite values of R. 



VII. RESULTS 



The results of our DMC calculations are presented in Figs. |4]-|10|. The error bars indicate 
our best estimate of the error due to statistical fluctuations as well as the contamination of 
higher energy states in the physical mass calculations. 

For the theory defined by C + we have calculated the vacuum expectation value and the 
physical mass over a range of coupling strengths[| from weak coupling (c ~ 0) to very strong 
coupling (c ~ 1). This calculation was done with J max = 10 and R = 3 which corresponds 
to (Jmax + I) 2 = 121 partial wave modes and a spherical surface of area 4nR 2 = 36ti. 
The vacuum expectation value is found to remain zero, i.e. the function g(xo,o) has a single 
maximum at xo,o = 0, throughout this range. Furthermore, as shown in Fig. |], the physical 
mass begins near unity at weak coupling and consistently grows larger as the coupling 
strength is increased to large values. We would expect the physical mass to vanish at the 
point where — > —(f) reflection symmetry is spontaneously broken. Our results clearly show 
that there is no phase transition in the theory corresponding to C + . 

On the other hand, our results do confirm the phase transition expected for the theory 
defined by £_. We calculated the vacuum expectation value and the physical mass using 



7 We work in units such that [i = 1 
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Jmax — 10 R— 3 



m 



♦ + + 



0.2 0.4 0.6 0.8 1 

c 

FIG. 4. m as a function of c for the theory defined by C + with J ma x = 10, i? = 3. 
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0.05 0.1 0.15 0.2 

c 

FIG. 8. m as a function of c for the theory defined by £_ with J m ax = 6, R = 2. 



m 




. 05 ^ ^ . 1 ^ .15 ^ 0.2 

c 

FIG. 9. m as a function of c for the theory defined by £_ with J m ax = 10, i? = 2. 
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0.8 



. 6 



m 



0.4 



0.2 




FIG. 10. m as a function of c for the theory defined by C_ with J m ax = 10, i? = 3. 



■/max = 6, R = 2; J max = 10, R = 2; and J ma x = 10, R = 3. The results are shown in 
Figs. |5|-|T0|. We find that at strong coupling the vacuum expectation value is zero but as we 
decrease the coupling strength the theory develops a nonzero vacuum expectation value and 
the physical mass vanishes. This clearly indicates that — > — reflection symmetry, while 
manifest at strong coupling, is spontaneously broken at weak coupling. 

The critical exponent f3 is defined by the behavior of the vacuum expectation value as we 
approach the phase transition from the broken-symmetry phase while the critical exponent 
v is defined by the behavior of the physical mass m as we approach the phase transition 
from the symmetric phase. Therefore, we have fit the results using the parametrized forms 



0" 



H 



<HC C rit - c 



and 



m 



K c ~ c crit) 1 



(60) 



The results of these fits are shown in Tables | and ||. 

The vacuum expectation value results are consistent with each other and with the recent 
calculations of (3 shown in Table |T|. Also note that c cr ^ < 0.3522. This is consistent with 
the finding of no phase transition for the C + theory as discussed below Eq. (|45|). 

The physical mass results are not as satisfying. We obtain values for v which are consis- 
tently higher than expected for the three-dimensional Ising universality class as seen by a 
comparison with the recent calculations of v shown in Table PTI} Furthermore, c CT \i is found 
to be lower than the value obtained from the vacuum expectation value calculations. These 
difficulties arise because of a systematic error in the calculation of the physical mass very 
near the phase transition due to finite-size effects. This systematic error is apparent in each 
of Figs. |8]-|TU| where the calculated values form a tail at small c. In the region of these tails 
the physical mass corresponds to a correlation length which is larger than the size of the 
system as given by R and so these points are of dubious validity. Therefore we must do our 
fits only for those points which are not too near the phase transition, i.e. for which m 
This restriction makes it problematic to extract an accurate value for v from fits to the 
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results. In the case of the vacuum expectation value calculations these finite-size effects are 
minimized as discussed below Eq. (|57). 



•Anax 


R 


a 


c crit 


P 


6 


2 


1.1(2) 


0.081(4) 


0.33(3) 


10 


2 


1.2(1) 


0.085(4) 


0.34(3) 


10 


3 


1.1(1) 


0.082(4) 


0.31(3) 




TABLE I. Fit parameters 


for vacuum expectation value calculations. 




<Anax 


R 


b 


c crit 


V 


6 


2 


3.2(4) 


0.064(6) 


0.72(7) 


10 


2 


3.8(7) 


0.077(7) 


0.69(9) 


10 


3 


3.7(6) 


0.074(7) 


0.78(9) 




TABLE II. Fit parameters for physical 


mass calculations. 




Ref. 


Method 





;/ 






MC 


0.3265(3)(1) 


0.6294(5)(5) 




§ 


HT 


0.32648(18) 


0.63002(23) 




I 


d = 3 exp. 


0.3258(14) 


0.6304(13) 




§ 


e-exp. 


0.3265(15) 


0.6305(25) 





TABLE III. Recent results for the critical exponents of the three-dimensional Ising universality 
class. 



VIII. SUMMARY 

We have used the formalism of spherical field theory to analyze the phase structure of 
three-dimensional 4 theory. This analysis is intended as a demonstration that this new 
nonperturbative method can be successfully extended beyond two-dimensional theories. We 
find that for fi 2 > there is no spontaneous breaking of the <ft ~^ ~ < f ) reflection symmetry 
while for fi 2 < this symmetry is broken at weak coupling and is restored as the coupling 
strength is increased. We have calculated the critical coupling and the critical exponents (3 
and v which characterize the phase transition in the latter case by analyzing the behavior of 
the vacuum expectation value and the physical mass as functions of the coupling. Our result 
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for (3, which was obtained from the vacuum expectation value analysis, is in agreement with 
recent calculations within the three-dimensional Ising universality class. However, our result 
for v is somewhat larger than the results of these calculations. This is most likely due to 
finite-size errors which affect the physical mass analysis. 

Several methods of improving this analysis come to mind. Our computations were per- 
formed with limited computing power: a 500 MHz Alpha workstation and a 350 MHz PC 
processor. Greater computing power could improve the analysis in two ways. First, the num- 
ber of iterations performed within the Monte Carlo algorithm could be increased thereby 
reducing the statistical errors in the computations and increasing the precision of the re- 
sults. Second, since the computation time scales roughly as J^ax a substantial increase 
in computing power would allow us to take larger values for J max and R in order to re- 
duce the systematic errors which arise due to these parameters being finite. Even with the 
limited computing power mentioned above, however, an improved analysis can be achieved 
with the method of periodic field theory jnj which uses a periodic-box mode expansion in 
place of the partial wave decomposition of spherical field theory and has the advantage of a 
time-independent hamiltonian. Improved calculations using this last approach are currently 
under way. 
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